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ABSTRACT 



We introduce a real-time decision support system 
which uses optimization methods, simulation, and 
judgement of the decision maker for operational 
assignment of military units to tasks and for tactical 
allocation of unit resources to task requirements. 

The system, named ARES for the Greek god of war, 
accommodates a high degree of detail in the logistics 
of unit movements during operations, yet separates the 
assignment and allocation activities in a fashion which 
naturally accommodates human intervention and 
j udgement . 

ARES is designed to assist the decision maker, not 
to replace him. ARES is demonstrated with a 
hypothetical scenario constructed for 14 Engineering 
Battalions of the Hellenic Army which are assigned 20 
tasks employing 25 resource types in repairing major 
damage to public works following a grate earthquake*. 

ARES is designed for use in real time, and quick 
data preparation is aided by the provision of standard 
task icons from published sources. 



* This hypothetical data was prepared prior to the 
earthquake in Kalamata near Athens on 13 September, 
1986, and exhibits uncanny, but coincidental 
resemblance to that real situation. 
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I. INTRODUCTION 



We introduce ARES, a prototypic system for real-time operational and 
tactical decision support. ARES is designed to quickly and effectively 
help respond to complex emergent problems in disaster relief, in the 
operational art and tactics of warfare, and in related multiperiod, 
large-scale employment of heterogeneous, substitutable resources 
restricted in availability and demand over time, over geography, and by 
organizational limitations . 

Although a great deal of work has been done in strategic modelling in 
many contexts, there is relatively little available modelling help beyond 
simple thumb rules for the time-pressed (operational or tactical) decision 
maker to translate strategic goals into logistical ly constrained 
operational and tactical plans (and the issues are different). The 
luxuries of hypothetical additional resources and the time to analyze 
their employment are just not available in the operational and tactical 
domains: operational and tactical decisions must be made quickly, and 

usually involve employing only resources actually available to perform 
whatever mission is at hand. 

The history of assignment and allocation models for planning emergency 
logistics extends back to some of the earliest work in linear programming, 
game theory, and their economic interpretation. We cite only a few of the 
references in this large body of literature. The seminal works by Dantzig 
and by Koopmans (both found in Koopmans 1951) are explicitly motivated by 
large-scale logistics problems. Karchere and Hoeber (1953) give early 
direction on the use of newly developed optimization technology in weapon 
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system planning and allocation, discussing substitutability of resources 
and choice of suitable objective functions. Geisler (1959) reports RAND’s 
first use of man-machine simulation of logistics support activities. 
Chaiken and Larson (1972) state some basic issues in logistic location and 
task assignment for emergency service vehicles: how many units should 

there be, where should they be located, who should they serve, and how can 
they be relocated to substitute for units not available. Kaplan (1973) 
redeploys divisible resources with linear programming. Fitzsimmons (1973) 
states a nonlinear response time model and uses pattern search to locate 
units well and allocate workload equitably. Swoveland, Uyeno, Vertinsky, 
and Vickson (1973) employ simulation and human interaction to set up a 
unit location problem as a quadratic assignment model which is solved with 
an elegant heuristic. Bracken and McGill (1974) formulate strategic force 
planning models as two-sided games solved with nonlinear programming. 
Bracken, Falk, and Karr (1975) apply multiperiod, two-person zero-sum 
games formulated to develop strategies for unit sortie allocations. 
Finally, Kolesar and Walker (1974) develop a multi-stage solution approach 
to unit and task assignment using set covering and transportation-1 ike 
integer linear programs which are used in real time by applying 
heuristics. 

Named for the Greek god of war, ARES is a proof prototype of a 
real-time decision support system. It employs optimization and simulation 
to capture and exploit a high degree of realism without demanding 
unreasonable amounts of data, or locking the decision maker out of the 
decision process. The intent is to provide quick credible advice with 
good global perspective at a cost no greater than the relatively myopic 
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decision methods now widely used. 

ARES accommodates enough detail to support realistic decisions, but 
not so much as to render the process useless. For the intended 
applications, the particular missions to be performed will not likely be 
known much in advance, but the generic types of missions are known and can 
be planned. ARES uses a taxonomy of prepared standardized icons for data 
describing possible missions. The idea is to help the decision maker 
quickly assemble a data scenario closely resembling the proximate 
situation from a menu of these standard icons. 

We characterize the mission at hand as a set of geographically 
dispersed tasks , each composed of partially-ordered sub- tasks requiring 
over time varying amounts of different resources . Organizational units , 
also geographically dispersed and each possessing a different endowment of 
resources, are to be assigned responsibility for the tasks. 

Responsibility for each task rests with only one unit at any given time. 

ARES consists of several models coordinated by a time interval 
simulator which also scales and manipulates scenario data in a fashion 
transparent to the decision maker. Two integer linear programs, a linear 
program, a georeference system, a mobility system, a decision-maker 
simulator, and extensive user interface and user override and control 
facilities complete the program suite. The models in ARES all use a 
standard data interface visible to the decision maker; this invites 
expansion with new models and features. 

A scenario is created from the attributes of available units and task 
attributes derived in large part from standard cataloged data icons for 
similar tasks. A georeference system is accomodated to generate distance 
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costs and delays in relocating and operating units. The decision maker 
may preview the scenario and modify data or manually pre-assign tasks and 
units as he sees fit. 

Operational assignment of tasks to units uses one of two integer 
programming models ( IPj^) or (IP). Good task aggregations for the unit 
assigned reduce unit relocation costs and match unit resource endowments 
with aggregated task resource requirements. Logistical considerations are 
paramount at this stage. 

The decision maker can review the operational assignments, modify them 
manually, or reject them outright and restate the conditions for the 
original operational assignment scenario. An acceptible set of 
operational assignments is passed forward to a tactical model. 

Tactical allocation of the resources of each unit to the requirements 
of its assigned tasks uses a linear programming model (GN). Substitutions 
among resources are permitted, although at reduced efficiencies in 
completing the tasks. Allocations recognize task priorities and the 
logistical effects of geographic proximity. In addition, unit efficiency 
in performing a particular task improves over time, and the sequence 
within tasks of resource requirements is considered. The result is a 
complete plan for each unit, showing what resources are to be used to 
fulfill each task requirement, and the efficiency with which operations 
are expected to be carried out. The allocation also determines which 
requirements will not be met in situations which overtax units. 

Finally, the decision maker is presented with a complete solution, 
which he can accept, or modify, or reject outright and repeat. Regardless 
of his action, ARES is designed to lend quick insight. The decision maker 
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uses his own judgement concerning non-quant i f ied factors, and he should 
gain a deeper understanding of the situation at hand from ARES. 
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II. OPERATIONAL ASSIGNMENT MODEL (IP) 



This integer program finds good aggregate assignments of tasks to 
units without explicit consideration of unit movement. 

Index Use 

i Tasks 

j Resources 

k Units 



Given Data 



d ik 



r . . , r . . 
-ij iJ 

a .. , a .. 
— jk jk 

“jk jk 

P i 



u . , u . 
—l 1 



Distance cost from unit k to task i 



f . . 
JJ 



ij 



Minimum , maximum resource j requirements of task i 



Minimum , maximum resource j employable by unit k 



h. 



Penalties for violating minimum , maximum resource limits 
Priority of task i (>0) 

Penalties for not assigning, double-assigning task i 
Substitution efficiency of resource j (>0) 

Sequence of resource j requirement within task i 
Consumption by task i of resource j from unit k 



i jk 

Decision Variables 

x^ Binary variable for assigning task i to unit k 

Formulation 

MIN 22 d.. x . . 

lk lk 
lk 



s . t . 2 
k 



x ik = 0- 1 ) : (u..u.) for all i (1) (GUB) 

o — — 

2 h. x.. = (a., ,a., ) ; (z ,z., ) for all j,k (2) 

ljk lk v - jk jk y v — jk jk y v ’ 



x ik - (0 ’ 1) 



for all i,k (3) 



(IP) 
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The notation = (r_ f r) ; (z^, z) indicates lower and upper ranges (r_, r) on 
row functional values with corresponding respective linear penalties per 
unit of violation (z,z); i.e., this is a goal program with linear 
penalties, an elastic integer program (Brown and Graves 1975). 

Constraints (1) encourage assignment of each task to some unit, and form a 
generalized upper bound (GUB) row set (Dantzig and Van Slyke 1967). 
Constraints (2) express the goodness of fit of task assignments with 
employable unit resources. Constraints (3) preclude fractional assignment 
of tasks to units. 

Consumption by task i of resource j from unit k is defined: 



-2 



h. .. 
ijk 



= r . . 

“ij 



-In f .. + (p.-l)/10 + d.. /a, + t.. 
jj J lk k lk 



(i) 



where is the speed of advance of a unit and t.^ is the number of 

periods that unit k has already been assigned task i. The rationale for 

the particular consumption function (1) amplifies the resource requirement 

r. . to account for the state of resource readiness f .., the task priority 
-iJ JJ 

p^ (making less important tasks appear more expensive), the logistic 
proximity of unit k and task i, d^/a^, and learning curve effect as a 
function of time since assignment, t.^. The data are scaled so that (1) 
is in conformity with policy guidance or the judgement of the decision 
maker. Alternate consumption functions may appeal in other situations. 

The distance costs d^ and penalties u . , u^ and z .^, z^ are expressed 
in commensurate units and deserve some thought by the modeller. For 
instance, z ^ may be interpreted as how much additional distance cost 
should be incurred before considering overtaxing maximum resource 
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employment for unit k; this is a direct expression of logistical 

efficiency. For simplicity in our tests, distance costs d^ are scaled by 
a policy parameter, z ^ and are part of the input script, u^ is 



defined as 100/p^ , and u^ equals 100. 
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III. OPERATIONAL ASSIGNMENT MODEL (IP L ) 



The purpose of this integer program is to find good movements of units 
to locations from which they will be assigned good aggregate groups of 
tasks to perform. 

Index Use 



i 


Tasks 


j 


Resources 


k 


Uni ts 


1 

Given Data 


Locations (assumed here to be collocated with tasks) 


d ik 


Distance cost from unit k to location 1 


s ii 


Distance cost from task i to location 1 


r jH 


Gross resource requirement j of task i performed 




from location 1 


r\s 

jlk 


Net resource availability j of unit k located at 1 



Decision Variables 



z ii 


Binary variable for assigning task i to location 1 


x lk 


Binary variable for moving unit k to location 1 
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Formulation 



MIN 22 g., z.. + 22 d.. x n 
. , to i 1 1 1 . , lk lk 



il 

2 

1 



kl 



z il 



2 

1 

2 

k 

-z . . +2 
11 k 



'? r j i 1 z i 1 * ^ a jlk x lk ' 

1 K 



= (l,l);(u ,u.) for all i (1) (GUB) 

for all k (2) (GUB) 

for all 1 (3) 

for all i , 1 (4) 

for all 1 , j (5) 



X lk = ( 1 • 1 ) * ( m « m ) 



X lk = (°> !) : ( m>m ) 



x lk = (0, 1) ; (m.m) 



Z il 



= { 0 . 1 } 

x i k = (0 - 1 > 



for al 1 i , 1 (6) 

for all 1 ,k (7) (IP L ) 



(IP^) uses the notation of (IP). Constraints (1) encourage assignment 
of each task to some location. Constraints (2) allow movement of each 
unit to some location. (A GUB row set is formed by constraints (1) and 
(2).) Constraints (3) attempt to restrict assignments so at most one unit 
is moved to any particular location, (4) require that a unit be moved to 
any location to which a task is assigned, and (5) attempt to match for 
each location and each resource an aggregate assignment of tasks which 
have gross resource requirements about equal to the net resource 
availability of the unit moved to that location to perform the tasks. 
Constraints (6) and (7) preclude fractional location of tasks and units. 

Gross resource requirement r^..^ represents the resource j estimated to 

be required at location 1 in order that task i actually receive . . 

^ J 



jil 



= r . . 
-ij 



-In f 



JJ 



(p.-l)/10 + 



E ii /0 , 



( 2 ) 



where a expresses the logistic radius of influence from any location; we 
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have used o = 100. The gross resource requirement (2) amplifies the 



resource requirement r.. in the same fashion as (1). 

lj 

Net resource availability a j^ represents the amount of resource j 

which unit k can deliver from its endowment a., forward to location 1. 

jk 

Unit k may be moving toward location 1 while supplying this net resource. 



a jlk = a jk f Jj ( “lk * C- a lk) (3) 

where is the fraction of time which the unit will spend at its 

destination location, and is the speed of advance. 

The distance costs d,, and g.,, and the penalties u. u., m, b and b 

IK 11 XI 

all render the same objective function units. In our work, m = 100, and 
u_ and u^ are defined as in (IP). The penalties for assigning too little 
(or too much) resource j to location 1 are b (or b) . We have used b = 0.1 
and b = 0.01 . 
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IV. TACTICAL ALLOCATION MODEL (GN k ) 



This linear program allocates resources to the tasks assigned to unit 
k. 

Index Use 

i Tasks 

j Resources 

w Work (resources required by tasks assigned to unit k) 

Given Data 



r . , r . 
-iw iw 



q. . q. 

— iw iw 



a , a .. 
“Jk jk 



Minimum , maximum work requirements w of assigned task i 



*Jk* z jk 



p i 
f . 

JW 



IW 

e . . 

1WJ 



Penalties for violating minimum, maximum work requirements 



Minimum , maximum resource j employable by unit k 



Penalties for violating minimum , maximum resource limits 
Priority of task i 

Substitution efficiency of resource j for work requirement w 

(> 0 ) 

Sequence of work requirement w in task i (>0) 

Efficiency of resource j used for work w on task i 
Decision Variables 



y. . Allocation of resource j to task i resulting in work w 
iwj J & 

Formulation 



MAX 222 e. . y. . 

. . . iwj J iwj 

St IWJ J J 

2 y i„j - <ilw’ 7 iw) : fora11 *•" (I) (GUB) 

J 

^ e iwj y i»j = (%k- a jk> ; %k- 7 jk> for aU j (2) 



y. . > 0 
iwj “ 



for all i,w,j (3) (GN k ) 
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(GN^) uses the notation of (IP). However, the dimensions of (GN^) 
discriminate between resources consumed, j, and work completed, w, 
explicitly representing substitutability of resources. Constraints (1) 
encourage allocation of sufficient work resources, while constraints (2) 
indicate the desired mix of employable unit resources. Constraints (3) 
require non-negative resource allocations. (GN^.) * s 311 e l ast i c 
generalized network (Brown and McBride 1984). 

Efficiency of resource j used for work w on task i is defined 



e . . 
iwj 



-In 

e 



f . 
jw 



+ (p . —1 )/10 + st /10 + d.,/a v 

X X W XK K 



( 4 ) 



where s? = max {0, s. - 
xw 1 xw 

allocation, and a. is the 
k 

employs the readiness and 
reduces efficiency if the 



t}, t is the last time period of this 

unit speed of advance. The efficiency (4) 

substitutability of resources via f . . st 

jw xw 

work w should not be started until period s. . 

xw 
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If model (IP) has been used for operational assignment. 



d ik = (°. d ik - CT k /2 }- ( 5 ) 

If unit k is to be advanced toward, or to, location 1 by model (IP^), 

d ik = max {0, d lk - o^/ 2} + g.j. (6) 

These distance costs d^ in (5) or (6), and penalties £L* W » Q ^ w . z_^, and 
z k are all intended to yield the same objective function units. For our 
tests, q . = 100/p . s + , and q. = 100. 

-Mw 1 1W 1W 
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V. CONSIDERATION OF LOGISTICS 



The efficiency with which a unit completes a task depends heavily upon 
logistical considerations. If a unit is remote from a task, or must be 
moved, its efficiency suffers. Figure 1 shows an idealized situation with 
unit k, task i, and location 1. 



r~ 



t 

I 

i 



“ 1 
i 






I 

L _ 



| 1 

I I 



I l 




unit 



Figure 1 • Idealized Geographic Logistical Scenario 



Model (IP) assigns tasks to units relying exclusively upon d^* (IP^) 

moves units to new locations and assigns tasks to be performed from these 
new unit locations. (IP^) reco 5 n ^ zes d^ and g.^. The distances d^ and 
dj^ are surrogates for logistical costs of assignment during the ensuing 
time period. Clearly, (IP^) is more appropriate for situations in which 
unit movements are expected, (IP) when they are not. (IP^) P rov id es the 
decision maker with a better opening gambit than does (IP) if the scenario 
involves significant initial redeployment of units. 

Tactical allocation models (GN) are given unit and task assignments 
and planned unit movements. Therefore, (GN) can allocate resources using 
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any logistic efficiency function of assigned distances, and of other 
attributes induced only from assignment such as weather effects, speed of 
unit movement, etc. (GN) can also substitute resources at somewhat 
reduced efficiency as well as prioritizing their immediate application. 
Given a fairly reasonable operational assignment, (GN) provides a 
high-resolution work plan with rich logistic detail and good face 
validity. 
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VI. AN EXAMPLE SCENARIO 



We demonstrate ARES with an example constructed for Engineering 
Battalions of the Hellenic Army. The mission scenario involves 20 tasks 
repairing major damage to public works following an earthquake. For our 
purposes, there are 14 units, each endowed with some of 25 resources. 
Figure 2 shows the units and tasks from the ARES input script. In the 
United States, the Department of the Army defines unit types in (1976), 
and task standards in (1973a). 



UNIT 


LABELS- LOCATIONS, AND PRIOR A 


SSIGNMENTS 






UU 


1 LABEL 


1 


1 X-COORD 


1 Y-COORD 


ISOA 


1 


1ST C0M8AT BN 


5-35 


20.20 


20.30 


150 


2 


2ND COMBAT BN 


5-35 


0A.50 


05.50 


150 


3 


3RD COMBAT BN 


5-155 


07.00 


17.25 


ISO 


A 


ATH C0M8AT BN 


5-155 


0 A . 90 


08.70 


150 


5 


1ST CONSTR BN 


5-115 


05.20 


1 A . 5 0 


100 


6 


2ND CONSTR BN 


5-115 


12.30 


18.75 


100 


7 


3RD CCNSTR BN 


5-115 


11.05 


15.60 


100 


B 


ATH CONSTR BN 


5-115 


C8.50 


15.25 


100 


9 


1ST AIRBCR BN 


5-195 


12.00 


05.00 


200 


10 


2ND AIR8CR BN 


5-195 


20.20 


20.20 


200 


I I 


1ST LIGH.EOUI .CO 


S-5B 


00. A0 


12. BS 


150 


12 


1ST ENG ARM BT 


5-1A5 


13.20 


08.60 


120 


13 


2ND ENG ARM BT 


5-1 AS 


09.25 


02 . 90 


120 


IA 


3RD ENG ARM BT 


S-IAS 


0B .50 


15.25 


120 


TASK 


LABELS. LOCATIONS, PRIORITIES, 


AND PRIOR ASS 


IGNMENTS 




TT 


1 LABEL 


1 


1 X-COORD 


1 Y-COORD 


IPRI 


1 


ADMIN BUILDING 


AA1051 


09.55 


05.90 


I 


2 


ADMIN BUILDING 


AA 1051 


09.60 


11.60 


I 


3 


ADMIN BUILDING 


AA1 101 


01.90 


06.95 


1 


A 


HOSPITAL 100 BED 


GH0I 11 


10.75 


07. A5 


1 


5 


HOSPITAL 200 BED 


GH02I I 


09.60 


11.60 


1 


6 


HOSPITAL 100 BED 


GH012I 


09.55 


05.90 


1 


7 


HOSPITAL 100 BED 


GH0 131 


01.90 


06 . 95 


1 


B 


RAILROAD BRIDGE 


B6I6A3 


09. 70 


05.90 


I 


9 


RAILROAD BRIDGE 


861512 


07.90 


09. A5 


I 


10 


ROAD BRIDGE 50* 


8SAI0I 


10.70 


07.20 


2 


11 


ROAD BRIDGE 100’ 


85A109 


09.70 


C5.B9 


1 


12 


ROAD BRIDGE 70' 


B5A10A 


08.20 


09.20 


2 


13 


ROAD 2.5 MILES 


853120 


10 .B0 


06 . 15 


2 


1A 


ROAD A. 7 MILES 


853122 


10.76 


07.50 


2 


15 


ROAD 5.5 MILES 


B53128 


0 9.60 


06.20 




16 


ROAD 6.8 MILES 


8 5 2 1 2 A 


02.25 


07.25 




17 


ROAD 6.0 MILES 


B 5 2 1 20 


10. B0 


06.95 


2 



CIPL) 
LL PP 



18 WATER T ANK-DI ST-SUP NOI 10.85 06.90 

19 WATER TANK-DIST-SUP N02 09.65 06.25 

20 WATER TANK-DIST-SUP NOS 09.10 10.85 



CIPJ(IPL) 



Figure 2 - Units and Tasks of Example 
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The georeference coordinates of units and tasks are given in Figure 1 
for the situation depicted in Figures 3, 4 and 5. 




0 10 ?0 



Figure 3: Initial Geographic Locations of Units. (Coordinates 

displayed are a georeference in common with the following 
f igures . ) 
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Figure 4: Earthquake Epicenter. 




Figure 5- Geographic Locations of Tasks. Geographic locations of 
damaged public works and earthquake epicenter are shown. 
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A georeference system is used to generate coordinate-to-coordinate 
distance costs, which appear in the ARES input script. 

The resource requirements for Task 1 ("ADMIN. BUILDING AA1051"), a 
disaster relief facility, are given in Figure 6 and the input script. 
Resource requirements such as these are available in standard engineering 
reference manuals for a wide variety of task types (for instance, see 
unclassified sources from the United States Department of the Army (1973a, 
1973b, and 1973c). We envision a taxonomy of standardized task data icons 
from which a particular set of requirements can be very quickly extracted 
and assembled for a scenario. The size of our resource requirements data 
base is modest but the resulting accuracy and level of detail are quite 
good. Better yet, data mobilization from a menu of such icons can be 
completed in minutes. 



RESOURCE 


LABELS AMD TASK 1 REQUIREMENTS 


MAN 


RR 


1 LABEL 1 


HOURS 


1 


ENGIN-P 10N-APREN-HLPER 


6648 


2 


SURVEYOR 


70 


z 


CARPENTER 


755 7 


4 


ELECTRICIAN 


940 


5 


PLUMBER 


1740 


6 


MASON 


UOO 


7 


STRUCTURE SPECIAL. 


0 


3 


HEAT-VENT1LAT SPECIAL. 


o 

o 

fj 


9 


WELDER 


0 


10 


PIPELINE 


0 


1 1 


CRANE-SHOVEL oper. 


rj 

o 

o 


12 


LOADER OPER. 


ISO 


13 


DOZER OPER . 


300 


14 


COMPRESSOR OPER. 


0 


IS 


DUMP TRUCK OPER. 


400 


16 


CONCRETE MACHINE OPER. 


0 


17 


GRADER OPER. 


329 


18 


CRUSHER O p ER . 


0 


19 


DITCH MACHINE OPER. 


300 


20 


ASPHALT SPECIAL. 


0 


21 


POWER ROLLER OPER. 


0 


22 


WATER DI3TRIBUT. OPER. 


0 


23 


POWER BOAT OPER. 


0 


24 


ROTARY TILLER OPER. 


0 


25 


SCRAPER OPER. 


0 



Figure 6: Resource Requirements of Task 1 
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The resources employable by Unit 1 ("1ST COMBAT BN"), a combat 
engineering battalion, are given in Figure 7 and in the input script. 
These resource endowments are in line with those given by the United 
States Department of the Army (1973a) with conversion to man hours from 
(1971). 



RESOURCE LABELS AMD UNIT I AVA 1 LA81 L I T I ES 



RR 


ILA8EL 1 


1 MIN 


IMAX 


IMIN PEN 


IMAX PEN 


z 


SURVEYOR 


A 05 


A50 


10 


10 


z 


CARPENTER 


1215 


1250 


10 


10 


A 


electrician 


203 


225 


10 


10 


5 


PLUMBER 


l A 18 


1575 


10 


10 


6 


3 

> 

o 

o 

z 


1215 


1350 


10 


10 


7 


STRUCTURE SPECIAL. 


0 


0 


10 


10 


8 


heat-ventilat special. 


0 


0 


10 


10 


9 


WELDER 


187 


208 


10 


10 


10 


PIPELINE 


0 


0 


10 


10 


11 


crane-shovel oper. 


1215 


1250 


10 


10 


12 


LOADER OPER. 


6165 


6850 


10 


10 


12 


DC2ER OPER. 


A 050 


A500 


10 


10 


1A 


COMPRESSOR OPER. 


1013 


1125 


10 


10 


15 


DUMR TRUCK OPER. 


10935 


12150 


10 


10 


16 


CONCRETE MACHINE OPER. 


203 


225 


10 


10 


17 


GRADER OPER. 


1620 


1800 


10 


10 


18 


CRUSHER OPER. 


0 


0 


10 


10 


19 


DITCH MACHINE OPER. 


0 


0 


10 


10 


20 


asphalt special. 


0 


0 


10 


10 


21 


POWER ROLLER OPER. 


0 


0 


10 


10 


22 


WATER D1STR1BUT. OPER. 


0 


0 


10 


10 


22 


POWER BOAT OPER. 


0 


0 


10 


10 


DA 


ROTARY TILLER OPER. 


0 


0 


10 


10 


25 


SCRAPER OPER. 


0 


0 


10 


10 




Figure 7: Resource 


Endowment of 


Uni t 


1 . 



The input script also includes for each task the sequence of resource 
requirements expressed as the first period when the resource is best 
applied, and for each resource its substitution efficiency for other 
resources . 

The scenario data constitutes about 1,000 records. However, these 
records derive from the unit, task, and resource definitions which are 
modest in number. 
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VII. DESIGN AND IMPLEMENTATION 



ARES is intended to help the decision maker, not to replace him. 
Figure 8 shows the functional structure of ARES. The design is biased 
toward interactive use with review and intervention options at each stage 
of operational assignment and tactical allocation. 
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INITIALIZE: 



NEXT_PERIOD: 
OP_ASSIGN : 



REVIEW_IP : 



TAC_ALLOC: 



UNIT-K: 



NEW_SCR I PT : 



RE V I E W_PER I OD : 



Define NEW_SCRIPT 

Redefine NEW_SCRIPT as OLD_SCRIPT 

Select Model (IP L ) or (IP) 

Read OLD_SCRIPT 

Generate and Solve (IP^) or (IP) 

Record task and unit assignments on ASSIGN_FILE 
Option to review assignments in ASSIGN_FILE 
either stop, 

or edit OLD_SCRIPT and GOTO OP_ASSIGN, 

or edit OLD_SCRIPT and/or ASSIGN_FILE and continue 

Read OLD_SCRIPT and store as SCRIPT 

Read ASSIGN_FILE and update SCRIPT assignments 

Select (GN^) , Generate and Solve 

Update SCRIPT resource requirements for work completed 
For next unit k REPEAT UNIT-K 

Update SCRIPT unit locations and distance costs 
Write SCRIPT as NEW_SCRIPT 
Option to review results 
either stop, 

or edit OLD_SCRIPT and/or ASSIGN_FILE 
and GOTO OP^ASSIGN 

or edit NEW_SCRIPT and GOTO NEXT_PERIOD 



Figure 8: ARES Functional Specification 
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ARES is implemented in FORTRAN H(Extended) and executes on an IBM 3033 
AP computer using the VM/CMS operating system. Input scripts are read 
from files which may be viewed and modified with a full-screen editor such 
as XEDIT. (Software copyrights IBM Corporation.) 

ARES uses the X-SYSTEM (Brown and Graves 1975) to solve (IP^), (IP) 
and (GN^) in real time. For each problem instance, problem generators 
directly convert input script data into an internal representation, the 
solver is invoked, and the solution is provided to a report writing 
program. ARES consists of a set of open subroutines and is executed with 
whatever preview, review, or other external interference is deemed 
desirable . 

We envision cyclic use and review at varying levels of detail as a 
mission progresses over time. Accordingly, input scripts include the 
beginning period and number of periods in the ensuing time interval, which 
intrinsically scales time-dependent input data to the desired level of 
aggregation. We have tested ARES manually and by replacing the decision 
maker with a simulation which performs " judgement review” of successive 
solutions over time. This permits totally automatic evaluation of 
complete mission scenarios, and avoids tedius manual effort in our 
research. (A single time interval may generate 15, or 20 thousand lines 
of solution detail at the scale of our example scenario.) 

The update of unit coordinate locations and distance costs is a simple 
surrogate for a more realistic and complicated georeference and mobility 
system. ARES estimates the direction and speed of advance of each unit 
during the time interval and relocates the unit. Then the distance costs 
are adjusted. If operating areas are known sufficiently in advance to 
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permit preparation of detailed georeference and mobility systems, ARES can 
accomodate the increased level of detail in real time (e.g. , Brown, Ellis, 
Graves, and Ronen 1987). The update can also be used to degrade, or to 
amplify unit resource endowments and effectiveness to modify task resource 
requirements, or to change any other data artifact, providing a rich 
modelling arena. 
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VIII. SCENARIO RESULTS 



ARES has been used in simulation mode to completely plan mission 
scenarios from start to finish. For the earthquake scenario. Figure 9 
shows the initial operational assignments of (IP^). 




0 10 20 

Figure 9: Initial Operational Assignments of Units. Directional 

vectors show the straight-line path and relative speed of 

advance a. . 

k 
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Figure 10 depicts the arrival of units to their initially assigned 



locations . 





Figure 10* Initial Operational Assignments of Units to Locations. 

Arrows show straight-line path of advance toward assigned 
locations . 

Without intervention by the decision maker, the scenario completed 
itself in 7 weekly intervals, requiring less than 2 minutes in a 1.2 
megabyte memory region. 

Face validity of the simulation solution is excellent. No manual 
intervention has been found to improve the solution. In fact, many manual 
attempts to coerce better assignments resulted in startling degradations. 
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The application of available resources, with allowable substitutions, 
is shown in Figure 11 for the 7 single-period time intervals required to 
complete the earthquake scenario. 



1 • t w • k gnd weak 3rd w»»k 4th woik 5th week 6th week 
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Figure 11: Resource Requirements and Work Completed. Each row 

represents a resource requirement over time-interval columns. 
The white bars depict resource requirements by time 
interval; the black bars show the relative fulfillment of the 
requirements. Broken bars are out of scale. From each time 
interval to the next the requirements are reduced by the work 
completed and amplified by new sequence-dependent 
requirements. In this scenario, 7 weekly time intervals are 
required to complete all tasks. 
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IX. COMPUTATIONAL EXPERIENCE 



Extensive computational experience reveals that the operational 
assignment models (IP) and especially (IP^) are most difficult to solve at 
the beginning of a scenario, and get progressively easy in later time 
intervals. The size of these models varies with the number of mandated 
assignments, impossible assignments, and the non-zero density of resource 
availabilities and remaining requirements. (IP) typically has about 340 
constraints, 268 binary variables, and 6,200 non-zero consumption 
coefficients. The linear program continuous relaxation can be generated 
and solved in about 5 seconds, and an optimal binary solution is achieved 
in another second, or so. 

(IP^) has about 1,000 constraints, 645 binary variables, and 8,000 
rather unwieldy non-zero gross resource requirement and net resource 
availability coefficients. 

The linear program continuous relaxation of ( IP^) proved impossible to 
solve by direct assault. Prior work by Brown and Graves for Bausch (1982) 
on large-scale set partitioning problems and later refinements by Brown, 
Graves and Ronen (1987) suggested an alternate means of at tack: a problem 

cascade. 

Briefly, the rows of constraints and columns of variables are 
lexicographically sorted to place short rows first accompanied by other 
rows and columns with intersecting non-zero coefficients, and longer rows 
later with their own intersecting rows and columns. 

The problem cascade proceeds by activating a set of constraints, 
relaxing all other constraints, and activating a set of variables, fixing 
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all other variables to their last-known values. This problem is solved, 
the new values of the active variables recorded, and another problem 
specified in the building problem cascade. The last problem in the 
cascade activates all constraints and variables (precisely the problem 
found intractable above) and solves i t by starting with an advanced 
solution recorded from the last-known values of variables solving previous 
problems in the cascade. 

(IP^) resisted even the problem cascade until a new heuristic cascade 
strategy was adopted which activates the shortest 1/2 of constraints and 
their associated variables, then the shortest 3/4, then 7/8, and so forth 
until the last constraint is added and the problem is solved. Remarkably, 
this approach has been absolutely reliable and robust, while most others 
fail or prove unruly. 

Generation and complete problem cascade solution of the continuous 
relaxation of (IP^) now requires about 10 seconds. 

An acceptable binary solution to (IP^) is achieved in another second, 
or two. 

We do not routinely seek optimal binary solutions to (IP^), which we 
refer to as ’’perfect misfits”. The gross resource requirements and net 
resource availabilities in (IP^) are rough logistic estimates, calibrated 
by actual field experience but ultimately just approximate target 
performance levels. For interesting operational assignments (i.e., early 
in the scenario) there are simply no feasible solutions; the goal is to 
guess where to send units so that they can peremptorily cope with the 
mission at hand with maximal effectiveness. Accordingly, we accept in 
practice binary solutions which may be as much as 25% greater than an 
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optimal lower bound in total value, including constraint violation 
penalties. Experimentally, we have determined at additional computational 
cost (as much as 5 minutes per trial) that these binary solutions are 
actually almost always within a few percent of the true optimum. 

A decision maker can help ARES with its operational assignments or 
completely specify a solution with manual assignment features. Our 
experience suggests that the decision maker can express some 
non-quantif iable guidance in this fashion, but can not hope to apply a 
remotely competitive global perspective. Manual competition with ARES 
reveals that model computation effort is amply justified by the quality of 
operational assignments achieved. The operational assignment models, 
especially (IP^), produce solutions no decision maker is likely to 
discover. Some of these solutions have yielded remarkable insights. The 
initial operational commitment of units is arduous and crucial to mission 
success. (IP^) * s wort ^ computational investment. 

By contrast, the tactical allocation models (GN) are easy to solve 
even in the cases where heroic substitution of resources are required. 

The size of each (GN^) varies with the number of tasks assigned to the 
unit, and the non-zero densities of resource availabilities, remaining 
requirements, and allowable substitutions. For our scenario, a typical 
instance of (GN^) has about 70 constraints and 1,190 variables, and is 
generated and solved in less than 0.1 second. Stress tests with 525 
constraints and 12,500 variables require less than a second. 
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X. DISCUSSION AND CONCLUSION 



The subtlety of operational assignment has surprised us, as has the 
ease of detailed tactical allocation. Operational assignments are 
delicate decisions, and the success of entire missions appears to be very 

sensitive to minute details precisely the considerations a hard-pressed 

decision maker would likely overlook in haste. 

Extensive mechanisms have been provided in ARES to encourage manual 
review and coercion of solutions. However, there have been very few cases 
in which such guidance improved solutions and many instances in which 
minor manual adjustments of operational assignments inflicted great 
disruption. For example, some operational assignments of (IP^) 
"cross-locate” units in the sense that a pair of units will each be 
collocated with a task assigned to the other. This superficial blemish 
can easily be masked by manual intervention or by automated solution 
editing. Surprisingly, the removal of cross-locations frequently 
increases the logistic cost of the solution: there is a very delicate 

balance of logistic support of task cohorts assigned to specialized units. 
Cross-location can actually make a great deal of sense in practice. 

Manual intervention can work well in cases inviting human judgement. 
For instance, nearly completed tasks or tasks which have been in progress 
for long intervals can enjoy efficiencies not apparent to our models. The 
decision maker can easily declare tasks completed when minor requirements 
remain, or when it is clear that the models are unduly influenced by a 
minor requirement. 
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Operational assignments can be restricted so that units are not moved 
from their initial new locations until the work in their logistic 
influence has been completed. Surprisingly, this restriction is rarely 
needed in practice, and in those cases in which multiple relocations are 
indicated great efficiencies accrue to the mission as a whole. We view 
this insight as a strong validation of the modelling philosophy underlying 
ARES. 

Fortuitous design decisions to separate operational assignment and 
tactical allocation models, to decompose time intervals, and to couple the 
resulting restricted components with simulation and human intervention 
options have yielded more than the intended benefits. Our original 
motives were to capture as much reality as possible while still rendering 
models capable of quick, responsive solution. 

The decomposed design also naturally accomodates features which are 
otherwise difficult to provide. For instance, partial orderings within 
tasks can be introduced. Also, discussions with Professor Wayne Hughes 
have suggested the technical feasibility of campaign analysis, two-sided 
gaming, and force-on-force applications of ARES. In these contexts, the 
coupling with simulation enhances our capabilities enormously. 

ARES was originally designed for use on an IBM-PC. There is no 
compelling reason not to use this microcomputer, but we encountered a few 
practical limitations. An arbitrary 640 kilobyte memory region limitation 
and crippling errors in the FORTRAN compilers available for the IBM-PC 
present artificial conversion costs which we are not willing to bear. 

When these unfortunate shortcomings of the IBM-PC are repaired, conversion 
might be reconsidered. Calibration tests project solution times on the 
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order of 2 minutes per time interval on IBM-PC/AT with a math co-processor 
and internal clock speed of 8 megahertz. 
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